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Abstract We review recent developments concerning the hadronic light-by-light scattering contribution to 
the anomalous magnetic moment of the muon. We first discuss why fully off-shell hadronic form factors should 
be used for the evaluation of this contribution to the g — 2. We then reevaluate the numerically dominant 
pion-exchange contribution in the framework of large- Nc QCD, using an off-shell pion- photon-photon form 
factor which fulfills all QCD short-distance constraints, in particular, a new short-distance constraint on the 
off-shell form factor at the external vertex in g — 2, which relates the form factor to the quark condensate 
magnetic susceptibility in QCD. Combined with available evaluations of the other contributions to hadronic 
light-by-light scattering this leads to the new result a L b y L ; had _ (H6±40) x 10 -11 , with a conservative error 
estimate in view of the many still unsolved problems. Some potential ways for further improvements are briefly 
discussed as well. For the electron we obtain the new estimate a L b y L ; had = (3.9 ± 1.3) x 10~ 14 . 
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1 Introduction 



The muon g—2 has served over many decades as an 
important test of the Standard Model (SM) . It is also 
sensitive to contributions from New Physics slightly 
above the electroweak scale. For several years now 
a discrepancy of about three standard deviations has 
existed between the SM prediction and the experi- 
mental value, see the recent reviews Refs. [l|-|4| on 
the muon g — 2. The main error in the theoretical 
SM prediction comes from hadronic contributions, i.e. 
hadronic vacuum polarization and hadronic light-by- 
light (had. LbyL) scattering. Whereas the hadronic 
vacuum polarization contribution can be related to 
the cross section e + e _ — > hadrons, no direct exper- 
imental information is available for had. LbyL scat- 
tering. One therefore has to rely on hadronic models 
to describe the strongly interacting, nonperturbative 
dynamics at the relevant scales from the muon mass 
up to about 2 GeV. This leads to large uncertainties, 
see Refs. @, 0, Q for recent reviews on had. LbyL 
scattering. 

The still valid picture of had. LbyL scattering as 
proposed some time back in Ref. is shown in Fig.Q] 





(a) 



(h) 





(c) (d) 

Fig. 1. The hadronic light-by-light scattering 
contribution to the muon g — 2. 

There are three classes of contributions to the rele- 
vant hadronic four-point function (WW) [Fig.[lja)] 
which can also be understood within an effective 
field theory approach to had. LbyL scattering: (1) a 
charged pion and Kaon loop [Fig. Hfb)], where the 
coupling to photons is dressed by some form fac- 
tor (p-meson exchange, e.g. via vector meson dom- 
inance (VMD)), (2) pseudoscalar exchange diagrams 
[Fig-HJc)] together with the exchanges of heavier res- 
onances (/ ,ai,...) and, finally, (3) the irreducible 
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part of the four-point function which was modeled 
in Ref. Q and later works @, Q by a constituent 
quark loop dressed with VMD-type form factors 
[Fig. Hid)]. The latter contribution can also be viewed 
as a short-distance complement of the employed low- 
energy hadronic models. According to quark-hadron 
duality, the (constituent) quark loop also models the 
contribution from the exchanges and loops of heav- 
ier resonances, like n',a' ,fQ,p,n,..., if they are not 
explicitly included in the other terms. 

One can try to reduce the model dependence 
and the corresponding uncertainties by relating the 
hadronic form factors at low energies to results from 
chiral perturbation theory and at high energies (short 
distances) to the operator product expansion. In this 
way, one connects the form factors to the underly- 
ing theory of QCD. This has been done in Refs. [8|- 
1 1 31 ] for the numerically dominant contribution from 
the exchange of light pseudoscalars ir°,r] and r]'. In 
Ref. 11| also important short-distance constraints on 
the axial- vector pole contribution have been imposed. 

2 On-shell versus off-shell form factors 

It was pointed out recently in Ref. Q, that one 
should use fully off-shell form factors for the evalua- 
tion of the LbyL scattering contribution. This seems 



to have been overlooked in the recent literature, in 
particular, in Refs. 0, E[_l3;12|. The on-shell form 
factors as used in Refs. |lol.ll2| actually violate four- 
momentum conservation at the external vertex, as 
observed already in Ref. 11 1. 

For illustration, we consider the contribution of 
the lightest intermediate state, the neutral pion. The 
key object which enters the diagram in Fig. [He) 
is the off-shell form factor J 7 7r o« 7 » 7 * ((q 1 + q 2 ) 2 , q 2 ,q 2 ) 
which can be defined via the QCD Green's function 

(VVP) @,iEl 



d 4 xd 4 ye^ x+ "^(0\T{j^x)j v (y)P 3 (0)}\0) 



(<?1 +<?2) 2 

x J : ^*^^((q 1 + q 2 ) 2 ,q\,ql) + (1) 



up to small mixing effects with the states r\ and rf and 

neglecting exchanges of heavier states like 7r ', 7r°", 

Here (x) is the light quark part of the electromag- 
netic current and P 3 (x) = (^ipi^fs^-ip^ {x). Note that 
for off-shell pions, instead of P 3 (x), we could use any 
other suitable interpolating field, like d^A 3 ^ (x) or even 
a fundamental pion field tt 3 (x). 

The corresponding contribution to the muon g — 2 
may be worked out with the result [Io| 



LbyL;7T° _ 



d 4 q 1 d 4 q 2 1 

(2tt)4 (2tt)4 q 2 1 qUq 1 +q2n(P + qi) 2 ~rnl}[(p~ q2 y-ml} 



q\-ml 



{{qi+q2) 2 ,q 2 ,q 2 ) J> 7 . 7 ((< ?1 +< ? 2) 2 ,(gi+<7 2 ) 2 ,0) 



0?i+<7 2 ) 2 -m 2 



T 2 (qi,q2-,P) 



(2) 



where the external photon has now zero four- 
momentum. See Ref. [lOj for the expressions for T t . 
Note that for general form factors a compact three- 
dimensional integral representation for a^^' 71 " has 
been derived in Ref. 131] - 



Instead of the representation in Eq. ([2]), Refs. [10, 
considered on-shell form factors which would yield 

the so called pion-pole contribution, e.g. for the term 

involving T 2 , one would write 

Jv<V 7 .(m 2 ,<2 2 ,g 2 ) x ^r0 7 . 7 (m 2 , (q 1 + q 2 ) 2 ,0). (3) 

Although pole dominance might be expected to give 
a reasonable approximation, it is not correct as it was 
used in those references, as stressed in Refs. 0, [Hj]. 
The point is that the form factor sitting at the 
external photon vertex in the pole approximation 



J r 7T o 1 , 1 (ml,{q 1 +q 2 ) 2 ,0) for (<?i+<7 2 ) 2 + ™ 2 violates 
four-momentum conservation, since the momentum 
of the external (soft) photon vanishes. The latter re- 
quires .7v>* 7*7(^1 +q 2 ) 2 ,(qi +g 2 ) 2 ,0). Ref. ljj] then 
proposed to use instead 



F^,^(ml,ql,ql) x J^o^rn 2 ,m 2 ,0) . (4) 

Note that putting the pion on-shell at the exter- 
nal vertex automatically leads to a constant form 
factor^nven by the Wess-Zumino-Witten (WZW) 
term 14|. However, this prescription does not yield 



the pion- exchange contribution with off-shell form 
factors, which we calculate with Eq. ([2]). 

Strictly speaking, the identification of the pion- 
exchange contribution is only possible, if the pion is 
on-shell. If one is off the mass shell of the exchanged 
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particle, it is not possible to separate different contri- 
butions to the g — 2, unless one uses some particular 
model where elementary pions can propagate. In this 
sense, only the pion-pole contribution with on-shell 
form factors can be defined, at least in principle, in a 
model-independent way. On the other hand, the pion- 
pole contribution is only a part of the full result, since 
in general, e.g. using some resonance Lagrangian, the 
form factors will enter the calculation with off-shell 
momenta. 

3 Pseudoscalar exchange contribution 

After the observation in Ref. Q that off-shell 
form factors should be used, the numerically dom- 
inant pion-exchange contribution was reanalyzed in 
detail in our paper [l3| . First we derived a new QCD 
short-distance constraint on the off-shell pion-photon- 
photon form factor ^" 7r o* 7 » 7 » from Eq. ([I} at the ex- 
ternal vertex in had. LbyL scattering. It arises in 
the limit when the space-time argument of one of the 
vector currents in (VVP) approaches the argument 
of the pseudoscalar density [15j ]. 

In the chiral limit, assuming octet symmetry and 
up to corrections of order a s , one then obtains the 
relation [la 



lim j;o« 7 , 7 ((Ag 1 ) 2 ,(A(7 1 ) 2 ,0) : 



(5) 



where F is the pion-decay constant in the chiral 
limit and x IS the quark condensate magnetic sus- 
ceptibility in QCD in the presence of a constant ex- 
ternal electromagnetic field, introduced in Ref. [lfj: 
{0\qa-^q\0) F = ee q x(ipip)o F^„, with e„ = 2/3 and 
e d = —1/3. Note that there is no falloff in Eq. ([5]) in 
this limit, unless x vanishes. 

Unfortunately there is no agreement in the litera- 
ture what the actual value of x should be. Note that 
X actually depends on the renormalization scale p. 
Most recent estimates yield values x(^ = 1 GeV) ~ 



3 GeV 2 [17|, although other approaches give a 



much larger absolute value of x{^ — 0.5 GeV) 
-9 GeV" 



16l . |18(. While the running with p can 



explain part of the difference, it seems likely that the 
different models used are not fully compatible. 

In Ref. we then reevaluated the pion- 

exchange contribution using an off-shell form fac- 
tor Jvo* 7 * 7 , ((g! -(- ffe) 2 ,? 2 ,*? 2 ) m the framework of 
large- N n Q CD. In the spirit of the minimal hadronic 
Ansatz [19J for Green's functions in large- N c QCD, 
such a form factor had already been constructed in 
Ref. 15]. It generalizes the usual VMD form fac- 



tor and contains the two lightest multiplets of vec- 
tor resonances, the p and the p' (lowest meson dom- 




inance (LMD) +V). In contrast to the VMD ansatz, 
the LMD+V form factor fulfills all the relevant short- 
distance constraints derived in Refs. 
eluding the new one from Eq. ([5]). In Ref. [1 
assumed that the LMD/LMD+V framework is self- 
consistent, therefore the estimate ^ LMD = — 2/M v = 
—3.3 GeV -2 was used (with a typical large- N c un- 
certainty of 30%), which is compatible with other es 
timates 



17] 



Other model parameters are fixed by normalizing 
the form factor to the pion decay amplitude n° — >• 77 
and by reproducing experimental data [2o| | for the on- 



shell form factor J r ^o 7 » 7 (m 2 , — Q 2 ,0), see Ref. [13[ for 



all the details. Recently, BABAR 21] has published 
new data for this form factor which does not show the 
characteristic falloff for large Euclidean momentum, 
lim Q 2^ 00 J"^o 7 . 7 (to 2 ,-(5 2 ,0) ~ 1/Q 2 [22J. Some im- 



plications of this new experimental result have been 
been discussed in Ref. |23j. As shown in Ref. 24], 
using the BABAR data to fit some of the LMD+V 
model parameters, does, however, not change the fi- 
nal result given below. 

Varying all the LMD+V model parameters in rea- 
sonable ranges and adding all uncertainties linearly to 
cover the full range of values obtained with the scan 
of parameters, we get the new estimate [3] 



LbyL;7r° 



: (72 + 12) x 10" 



(6) 



As far as the contribution to et M from the ex- 
changes of the other light pseudoscalars r\ and 77' 
is concerned, a simplified approach was adopted in 
Ref. as was done earlier in other works 
We took a simple VMD form factor, normalized to 
the experimental decay width F(PS — > 77). In this 
way one obtains the results a^ byL;7) = 14.5 x 10 -11 and 
a LbyL;r,' = i2.5xl0^ 11 . Adding up the contributions 
from all the light pseudoscalar exchanges, we obtain 
the new estimate 11311 



LbyL;PS _ 



(99 + 16) x 10" 



(7) 



where we have assumed a 16% error, as inferred above 
for the pion-exchange contribution. 

For comparison, we have listed in Table [T] some 
evaluations of the pion- and pseudoscalar-exchange 
contribution to had. LbyL scattering by various 
groups. The model used by each group has also 
been indicated in the first column of the table, see 
the corresponding references for all the details (we 
used the abbreviations: ENJL = Extended Nambu- 
Jona-Lasinio model; HLS = Hidden Local Symmetry 
model; xQM = chiral quark model; FF = form factor; 
h 2 is one of the LMD+V model parameters). 
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Table 1. Results for the n°,r/ and if exchange contributions obtained by various groups. 



Model for .7^0-. q^ byL;7r X 10 11 q^ byL;PS X 10 11 



Point coupling 


+oo 


+oo 


modified ENJL (off-shell) [BPP] [8] 


59( 9 ) 


85(13) 


VMD / HLS (off-shell) [HKS, HK] [9] 


57( 4 ) 


83( 6 ) 


nonlocal xQM (off-shell) [DB] [25] 


65( 2 ) 




AdS/QCD (off-shell ?) [HoK] [26] 


69 


107 


LMD+V (on-shell, h 2 = 0) [KN] [10] 


58(10) 


83(12) 


LMD+V (on-shell, h 2 = -!0 GeV 2 ) [KN] [10] 


63(10) 


88(12) 


LMD+V (on-shell, constant FF at external vertex) [MV] fll] 


77( 7 ) 


114(10) 


LMD+V (off-shell) [N] \M 


72(12) 


99(16) 



Our results for the pion and the sum of all pseu- 
doscalar exchanges are about 20% larger than the val- 
ues in Refs. 0, Q which used other hadronic models 
that presumably do not obey the new short-distance 
constraint from Eq. ([5]) and thus have a stronger 
damping at large momentum. Within the non-local 
xQM used in Ref. |25| there is a strong, exponential 
suppression for large pion virtualities. According to 
Ref. [itij . the estimate with the AdS/QCD model has 
an error of at most 30%. On the other hand, our re- 
sult is smaller than the pion- and pseudoscalar-pole 
contribution calculated in Ref. 11 1. Since only the 



pion-pole contribution is considered in Ref. , their 
short-distance constraint cannot be directly applied 
to our approach. However, our ansatz for the pion- 
exchange contribution agrees qualitatively with the 



short-distance behavior of the quark-loop derived in 

BEL 



Ref. 



11] 



Note, 



see the discussion in Refs. 
however, that the numerical value for the pion-pole 
contribution listed as [MV ] in Table Q] should rather 
be 80 x 10- 11 , see Refs. 0,E!,[25|. 



4 Summary of other contributions 

In Table [2] we have collected the results for all 
the contributions to had. LbyL scattering according 
to Fig. [1] obtained by various groups in recent times, 
including some "guesstimates" for the total value. In 
the following, we highlight the main features of the 
numbers given and point out some critical issues re- 
garding each contribution. A more detailed discus- 
sion can be found in Ref. Q. 



Table 2. Summary of the most recent results for the various contributions to aJ; byL;had xl0 11 . The last column 
is our estimate based on our new evaluation for the pseudoscalars and some of the other results. 



Contribution 


BPP [8] 


HKS, HK [9] 


KN [10] 


MV [11] 


BP [5J, MdRR [1] 


PdRV [6] 


N [13], JN [3] 


n°,r),r)' 


85+13 


82.7+6.4 


83+12 


114+10 




114+13 


99+16 


axial vectors 


2.5+1.0 


1.7+1.7 




22+5 




15+10 


22+5 


scalars 


-6.8+2.0 










-7+7 


-7+2 


n,K loops 
ir,K loops 
+subl. N c 
quark loops 


-19+13 
21+3 


-4.5+8.1 
9.7+11.1 




0+10 




-19+19 
2.3 


-19+13 
21+3 


Total 


83+32 


89.6+15.4 


80+40 


136+25 


110+40 


105+26 


116+39 



As one can see from Table [21 the different mod- 
els used by various groups lead to slightly different 
results for the individual contributions. The final re- 
sult is dominated by the pseudoscalar exchange con- 
tribution, which is leading in large- Nc, but sublead- 
ing in the chiral counting. The other contributions 
are smaller, but not negligible. Furthermore, they 
cancel out to some extent, in particular the dressed 
pion and Kaon loops and the dressed quark loops. 

In Ref. [111, new QCD short-distance constraints 
were derived for the axial-vector pole contribution 



with on-shell form factors J- Al * 1 * at both vertices. 
A huge enhancement of a factor of ten was observed 
compared to the earlier estimates in Refs. which 
assumed nonet symmetry for the states Oi,/i and f[. 
It was shown that the result is very sensitive to the 
mass of the exchanged axial-vector resonance. Since 
the form factors include light vector mesons like the p, 
this leads to a smaller effective mass of the exchanged 
resonance, compared to M A ~ 1300 MeV. The result 
is also sensitive to the mixing of the states /i and f[. 
The result given in Table [2] corresponds to ideal mix- 
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ing. If /i is a pure octet state and f[ a pure singlet, 



the final result goes down to a 



LbyL;ai,/i./{ _ 



17X10" 11 . 



11 



is an important 



The procedure adopted in Ref. 

improvement over Refs. 0, Q and we have therefore 
taken the result for the axial-vectors from that ref- 
erence for our final estimate for the full had. LbyL 
scattering contribution. This despite the fact that 
only on-shell form factors have been used in Ref. 11 1 . 
As we argued above, we think that one should use 
consistently off-shell form factors at the internal and 
the external vertex. 

Within the ENJL model used in Ref. [1] , the scalar 
exchange contribution is related via Ward identities 
to the constituent quark loop. In fact, Ref. |9| ar- 
gued that the effect of the exchange of scalar reso- 
nances below several hundred McV might already be 
included in the sum of the (dressed) quark loops and 
the (dressed) pion and Kaon loops. Such a potential 
double-counting is definitely an issue for the broad 
sigma meson / (600). It is also not clear which scalar 
resonances are described by the ENJL model used in 
Ref. [H]. The parameters were determined from a fit 
to various low-energy observables and resonance pa- 
rameters, among them a scalar multiplet with mass 
M s = 983 MeV. However, with those fitted parame- 
ters, the ENJL model actually predicts a rather low 
mass of M| NJL = 620 MeV. 

The (dressed) charged pion- and Kaon-loops from 
Fig- IUb) yield the leading contribution in the chi- 
ral counting, but are subleading in N c - We note 
that the result without dressing (scalar QED) is ac- 
tually finite: a^ byL;7r± = -46 x 1CT 11 . The dressing 
with form factors then leads to a rather huge and 
very model dependent suppression (compare the re- 
sults for Refs. 1) and Q in Table [2]), so that the final 
result is much smaller than the one obtained for the 
pseudoscalars. This effect was studied in Ref. 11 1 
for the HLS model used in Ref. [9], in an expansion 



in (m^/Mp) 2 . Ref. [ll| observed a large cancellation 
between the first few terms in the series and the ex- 
pansion converges only very slowly. The main reason 
is that typical momenta in the loop integral are of 
order /x = 4m T « 550 MeV and the effective expan- 
sion parameter is fi/M p . The authors of Ref. 11 1 
took this as an indication that the final result is very 
likely suppressed, but also very model dependent and 



that the chiral expansion looses its predictive power. 
The pion and Kaon loops contribution is then only 
one among many potential contributions of 0(1) in 
N c and they lump all of these into the guesstimate 
a LbyL,Ar c _ _j_ j^q-j x io -11 . However, since this es- 
timate does not even cover the results for the pion 
and Kaon loops given in Refs. HE], we think this 
procedure is not very appropriate. 

The (dressed) constituent quark loops from 
Fig. Old) are also leading in large- N c - The result with 
point-like couplings is finite: a LbyL ;qU ark S = 62 x 10* 11 . 
The dressing with form factors then leads again to a 
large and very model dependent suppression of the 
final result, compare Refs. @] and @] in Tabled 

In the recent review Q the central values of some 
of the individual contributions to had. LbyL scatter- 
ing were adjusted and some errors were enlarged to 
cover the results obtained by various groups which 
used different models, see Table [2] Finally, the er- 
rors were added in quadrature. Maybe the result- 
ing small error masks some of the uncertainties we 
still face in had. LbyL scattering. Note that the 
dressed light quark loops are not included as a sep- 
arate contribution in Ref. Q (only the contribution 
from a bare c-quark is included in Table [2]). The light 
quark loops are assumed to be already covered by us- 



ing the short-distance constraint from Ref. ll| on 
the pseudoscalar-po/e contribution. Although numer- 
ically the final estimate from Ref. 



11 



is very close 

to our result given in Table [2j in view of the inter- 
pretation given for this term in the Introduction, we 
do not see any reason, why the contribution from the 
dressed quark loops should be discarded completely. 
At least in large- N c QCD, only the sum of all reso- 
nance exchanges should be dual to the quark loops. 

5 Conclusions 

Combining our result for the pseudoscalars with 
the evaluation of the axial-vector contribution in 
Ref. 11 1 and the results from Ref. Q for the other 
contributions, we obtain the new estimate 0, [l3| 



a LbyL;had = ( - 116±40 ) xl0 - 



(8) 



for the total had. LbyL scattering contribution to the 
anomalous magnetic moment of the muonQ The vari- 



LbyL;7r° 



(2.98 ±0.34) X 10 -14 0. Note that the naive rescal- 



* Applying the same procedure to the electron, we get a e 

(me/nifj,) 2 a^ byL,7r = 1.7 X 10~ 14 yields a value which is almost a factor of 2 too small. Our cs- 

LbyL;,? _ „ 1n -14 ,„j LbyLjr;' 



ing Oj l>yL, ' r (rescaled) 



timatcs for the other pseudoscalars contributions are a^ uyLJ ''' = 0.49 X 10 — 14 and a e y ,r> = 0.39 X 10 -14 . Therefore we get 
^LbyL.PS _ (3 g-|-o.5) x 10 -14 . Assuming that the pseudoscalar contribution yields the bulk of the result of the total had. LbyL 
scattering correction, we obtain a LbyL.had _ ^ g± \ 3) x 10 -14 , with a conservative error of about 30%, see Ref. This value 
was later confirmed in the published version of Ref. [(J where a leading logs estimate yielded a LbyL ; had _ ^ 5 -j- ^ .0) X 10 — 14 . 
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ation of the results for the individual contributions 
listed in Table [5] reflects our inherent ignorance of 
strong interaction physics in had. LbyL scattering. 
One can take the differences between those values as 
an indication of the model uncertainty and, to be con- 
servative, all the errors have been added linearly, as 
was done earlier in Refs. [l|, 0, lil. Eo||. 

Certainly, more work on the had. LbyL scattering 
contribution is needed to fully control all the uncer- 
tainties, in particular, if we want to fully profit from 
a potential future g — 2 experiment with an expected 
error of about 15 x 10 -11 27f . Maybe at some point 



we will get an estimate from lattice QCD [28j . al- 
though the relevant QCD Green's function (WW), 
to be integrated over the phase space of three off-shell 
photons, is a very complicated object. 

In the meantime we suggest the following way 
forward [3j. It is very important to have a unified 
framework (hadronic model) which deals with all the 
contributions to had. LbyL scattering. A purely phe- 
nomcnological approach would be to use some reso- 
nance Lagrangian where all couplings are fixed from 
experiment. Since such Lagrangians are in general 
non-renormalizable it is, however, not clear how to 
achieve a proper matching with QCD at short dis- 
tances. Such a matching can be achieved within 
the large- N c framework, however, the corresponding 
resonance Lagrangians in general contain many un- 
known coefficients and it will be difficult to fix all 
of them theoretically or experimentally. In any case, 
in both of these approaches any additional experi- 
mental information on various hadronic form factors 
would be very useful to constrain the theoretical mod- 
els. In this respect, e + e~ colliders running at energies 
around 0.5—2 GeV could help to measure some of the 
form factors relevant for had. LbyL scattering [29| . 

/ would like to thank the organizers of PHIPSI09 
for the invitation and their financial support. I am 
grateful to F. Jegerlehner for many helpful discus- 
sions and numerous correspondences. This work was 
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